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ABSTRACT 

WASP- 13b is a sub-Jupiter mass exoplanet orbiting a G1V type star with a period of 4.35 days. 
The current uncertainty in its impact parameter (0 < b < 0.46) resulted in poorly defined 
stellar and planetary radii. To better constrain the impact parameter we have obtained high 
precision transit observations with the RISE instrument mounted on 2.0 m Liverpool Tele- 
scope. We present four new transits which are fitted with an MCMC routine to derive accurate 
system parameters. We found an orbital inclination of 85.2 ± 0.3 degrees resulting in stellar 
and planetary radii of 1.56 ± 0.04 R Q and 1.39 ± 0.05 Rj up , respectively. This suggests that 
the host star has evolved off the main-sequence and is in the shell hydrogen-burning phase. 
We also discuss how the limb darkening affects the derived system parameters. With a den- 
sity of 0.17 pj, WASP- 13b joins the group of low density planets whose radii are too large 
to be explained by standard irradiation models. We derive a new ephemeris for the system, 
T = 2455575.5136 ± 0.0016 (HJD) and P = 4.353011 ± 0.000013 days. The planet equi- 
librium temperature (T equ = 1500 K) and the bright host star (V = 10.4 mag) make it a good 
candidate for follow-up atmospheric studies. 

Key words: methods: data analysis - methods: observational - stars: planetary systems - 
stars: individual (WASP- 13) -techniques: photometric 
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The disc overy of the first short-period transiting planet HD 
209458b JCharbonneau et alj |2000| : iHenrv et alj 120001) opened a 
new field of exoplanetary research. Their advantageous geome- 
try allow us to estimate accurate planetary mass and radii. The 
bulk density of the planet provide s information on its composition 
(Guillot 2005; Fortney et al. 2007) and places constraints on plan- 
etary structure and formation models. 

Interestingly, many of the short- period planets are b loated 
relative to baselin e models ( e.g. Bodenh eimer et all 1200 3l : 
iFortnev et all [2007), which assume a hydrogen/helium planet 
contracting under strong stellar irradiation. We define the ra- 
dius anomaly of exoplanets as the difference between the mea- 
sured radius and theoretical radius. Negative radius anomalies, i.e. 
denser planets, can be explaine d by different planet composition 
or/and the presence of a core dGuillot et alll2006l : I Burrows et alj 
2007). However, the majority of the hot giant planets have low 
density (i.e. radius anomalies larger than zero) that are harder 
to reconcile. Several theories have been presented to explain 
the positive radius anomalies of hot exoplanets: these include 
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tidal h eatin g I B oden heimer et all 1 200 lh ■ kinetic heating due to 
wind s dGuillot & ShowmanT 20021) , enhanced atmo spheric opac- 



ities iBurrows et alj l2007t) and ohmic dissipation iBatygin et alj 
l201lh . Although these theories can ex plain individual sys tems, 
only the recent ohmic dissipation theory dBatygin et al ] l201lh . can 
succ essfully explain all of the known exoplanets radius anoma- 
lies dLaughlin et alj|2oTlh . The lowest density hot Jupiters, such as 
WASP- 13b, are particularly important to shed light on the physical 
mechanisms that lead to inflated planets. 

For bright transiting systems, further insight into their phys- 
ical properties can be obtained through followup observations. 
Transmission spectroscopy, which consists of measuring the stellar 
light filtered through the planet's atmosphere during transit, pro- 
vides information into exoplan et atmospheres 1 Charb onneau et all 
2002; [Vidal -Madi aret al.l2003f) . Observation of secondary eclipses 
(i.e. occultations) offers the potential for directly measuring plan- 
etary emission spectra (e.g. Deming et al] 1 20051 : iGrillmair et alj 
2008; ICharbonneau et all l2008h . Finally, the ability to measure 
the sky-projected angle between the stellar rotation axis and 
planetary orbit, through the Rossiter-McLaughlin effect jRossiteJ 
1 19241 ; [McLaughlin 19 2^), can yield clues about the formation and 
migrati on proces ses ( Fabr vckv & Wrrinl 120091 ; iTriaud et al] l20ld : 
IWinn et alJliciToh . 
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Followup opportunities benefit from accurate planetary and 
stellar parameters. However, obtaining high signal-to-noise tran- 
sit observations from the ground is difficult. Consequently even 
some of the planets with bright host stars are lacking good quality 
light curves, and hence have poorly de termined planetary param- 
eters. This is the case for WASP-13b dSkillen et ai]|2009l) where 
it was not possible to accurately constrain the impact parameter 
(0 < b < 0.46) due to the lack of a high precision transit light 
curve. 

WASP-13b is a low-mass p lanet with M v = 0.46 ± 0.06 
Mj up in a 4.3 day circular orbit JSkillen et alj|2009l) . Its host star 
is G1V type with M» = 1.03 ± 0.10 M , T eff = 5826 ± 100 K, 
logg = 4.04 and solar metallicity. WASP- 1 3b was discovered i n 
2006-2007 by the SuperWASP-north survey jPollacco etaT]|2006h . 
As mentioned above, the impact parameter is not well constrained 
and two possible sets of system parameters are presented in the 
discovery paper for two possible impact parameters; 6 = and 
b = 0.46. This results in large uncertainties in both stellar (R, ~ 
1.2 - 1.34R ), and planetary radii (Rp ~ 1.06 - 1.21). The 
adopted solution b — 0.46 suggests that the host star has evolved 
off the main-sequence to the shell hydrogen burning phase. There- 
fore, this system is interesting in the context of exoplanet evolution, 
and with a bright host star (V=10.4 mag) is a good candidate for 
followup studies. 

In this paper, we present four high precision transit observa- 
tions of WASP- 13b which are described in Section 2. Our data al- 
low us to determine the inclination of the system and directly de- 
rive the mean stellar density and hence, better constraint the stellar 
and planetary radii. We discuss our transit model in Section 3 and 
present the updated parameters of the system in Section 4. In Sec- 
tion 5 we discuss the discrepancy between theoretical and empirical 
limb darkening coefficients. Finally, we discuss our results in Sec- 
tion 6. 



2 OBSERVATIONS 

WASP- 13b was o bserved wi t h RI SE dSteele et alj 120081 ; 
iGibson et al] 120081 ; iBarros et all 1201 lh mounted on the 2.0m 
Liverpool Telescope on La Palma, Canary Islands. RISE is a 
frame transfer e2v CCD with a pixel scale of 0.54 arcsec/pixel 
that results in a 9.4 x 9.4 arcmin field-of-view. It has a wideband 
filter covering ~ 500-700 nm which corresponds approximately to 
V+R. For exposures longer than 1 second, RISE has a deadtime of 
only 35 ms. To reduce the deadtime to the minimum (Barr os et al.l 
1201 lh we used an exposure time of 1 second which required a 
defocussing of -1.0 mm to avoid saturation of the bright com- 
parison star (8.7 mag). Defocussing i s also important to decrease 
systematic noise due to poor guiding l lBarros et alj|201 lh . 

The data wer e reduced using the ULTRACAM pipeline 
dDhillonetaT1l2007h which is optimized for time-series photom- 
etry. Each frame was bias subtracted and flat field corrected. We 
performed differential photometry relative to an ensemble of com- 
parison stars in the field, confirmed to be non-variable. For each 
observation we sampled different aperture radii and chose the aper- 
ture radius that minimised the noise. The photometric errors in- 
clude shot noise, readout and background noise. 

On 2009-01-29 we obtained an ingress of WASP- 13b and 
an egress on 2009-05-05. We obtained full transits of WASP-13b 
on 2010-02-03 and 201 1-01-13. Unfortunately, both observations 
were degraded by some clouds. Observations with transparency 
lower than 94% were clipped. On 2010-02-03 the transit obser- 
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Table 1. WASP- 13 observation log. 



Date 


number exposures 


FWHM 


aperture radius 






pixel 


pixel 


2009-01-29 


4400 


17 


23 


2009-05-05 


10780 


21 


24 


2010-02-03 


17000 


17 


22 


2011-01-13 


16000 


24 


27 



vation was interrupted after ingress due to the derotator reaching 
its limiting position and the two parts of this transit are treated as 
separate observations. This observation shows systematics near the 
rotator interruption and is affected by clouds before the egress. We 
experimented removing it from the final analysis and concluded 
that its inclusion does not bias the results. On the last 20 minutes 
of the 2011-01-13 observation the position of the star in the CCD 
varied by 10 pixels in the x and 15 pixels in y which caused an addi- 
tional trend in the light curve. Due to its effect on the normalisation 
of the light curve, the last 20 minutes of the 201 1-01-13 observation 
were clipped. The number of exposures taken, estimated FWHM 
and aperture radius used for each observation are given in Table [T] 
We also included in our analysis, the previously published transit 
of WASP-13 taken with 0.95 m James Gregory Telescope (JGT) at 
St Andrews University on 2008-02-16 ( Skil len et alj2009h . 

The final high precision transit light curves for WASP- 13b are 
shown in Figure [T] along with the best-fit model described in Sec- 
tion |3.2| We overplot the model residuals and the estimated uncer- 
tainties, which are discussed in Section [3TI To illustrate the quality 
of our results we also show the phase folded weighted combina- 
tion of all the light curves in Figure [2] Around -0.01 in phase the 
only high precision data is from the 2011-01-13 light curve that 
shows extra scatter, is below the transit model. The remaining tran- 
sit phases are better sampled and hence have smaller errors. 



3 DATA ANALYSIS 
3.1 Photometric errors 

As mentioned above the initial photometric errors include only shot 
noise, readout and background noises, which underestimates the 
true errors. To obtain a more reliable estimate of the errors we be- 
gin by scaling the errors of each light curve so that the reduced 
X 2 of the best fitting model is 1.0. This resulted in the multipli- 
cation of the errors by 1.01, 1.55, 2.46, 2.68, 2.26 and 2.79, for 
the 2008, Jan 2009, May 2009, ingress 2010, egress 2010 and 201 1 
light curves, respectively. For the RISE observations, this indicates 
that the photometric errors are higher than expected from photon 
noise, probably due to the non-photometric conditions of these ob- 
servations. 

Exoplanet transit observations are also affected by time cor- 
related noise which can lead to underestimate d parameter uncer- 
tainties jPontetall 2006: Cart er & Winr]|2009h . The residuals dis- 
played in FigureQjshow that the photometric noise is dominated by 
white noise. Nevertheless, we included time-c orrelated noise fol- 
lowing the procedure from Gillon et al.l d2009l) . Using the residu- 
als of the best fit model, we estimated the amplitude of the red 
noise, ay, to be 100 ppm (Jan 2009), 50 ppm (May 2009), 250 ppm 
(ingress 2010), 150 ppm (egress 2010), and 200 ppm (2011). For 
the 2008 light curve, the red noise was found to be negligible. These 
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Figure 1. Transit observations of WASP- 13b. From top to bottom in chrono- 
logical order; JGT 2008 February 16, RISE 2009 January 01, RISE 2009 
May 05, RISE 2010 February 03 and RISE 2011 January 13. We super- 
impose the best-fit transit model and also show the residuals for each light 
curve at the bottom of the figure. The RISE data are binned into 30 second 
periods, and bins displaced vertically for clarify. The individual RISE light 
curves plotted here are available in electronic form at CDS. 



errors were added in quadrature to the rescaled photometric errors 
and were used in the final Markov Chain Monte Carlo (MCMC) 
chains. 



3.2 Determination of system parameters 

To determine the planetary and orbital parameters, we fitted the 
four RISE light curves and the original JGT light curve of WASP- 
13b simultaneousl y. We used the same proc edure as for WASP-21b 
( lBarrosetai]|201lh which is based on the lMandel & Agoll ( 2002) 
transit model parametrised by the normalised separation of the 
planet, a/7?*, ratio of planet radius to star radius, Rp/R*, orbital 
inclination, i, and the transit epoch, To, of each light curve. We ac- 
count for a linear trend with time for each light curve parametrised 
by the out of transit flux and flux gradient. Due to the short or- 
bital period and age of the system we expect the orbit has circu- 
larised (see below), hence a circular orbit was assumed. The transit 
model is coupled with an MCMC proced ure to obtain a ccurat e pa- 
rameters and uncertainties as outlined in iBarros etalJfeOllh . We 
included the quadratic limb darkening coefficients (LDC) from the 
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Figure 2. Phase-folded combined transit light curve of WASP- 13b. The data 
was optimally binned in 150 phase bins for clarity. We superimpose the 
best-fit transit model. Note that here we do not account for the different 
limb darkening of the RISE and the JGT light curves. 



models of iHowarthl fcOl 1 j) forT e// = 6000 K, logg = 4.0 and 
\M/ H] = 0.0 whi c h wer e the closest tabulations to the parameters 
from lSkillenetalJ d2009l) : for the RISE filter V+R, a = 0.4402 
and 6 = 0.2394 and for the R filter of JGT, a = 0.4121 and 
b = 0.2312. We initially kept the limb darkening parameters fixed 
during the fit. We treated the two parts of the 2010-02-03 observa- 
tions as separate transits to correctly account for the normalisation 
but restricted them to have the same To. Including the two linear 
normalisation parameters for each light curve, a total of 20 param- 
eters were fitted. Besides the linear normalization, no extra trends 
were removed from the light curves. 

We computed seven MCMC chains, each of 500 000 points 
and different initial parameters. The initial 20% of each chain that 
corresponded to the burn in phase were discarded and the remaining 
parts merged into a master chain. We define the best fit parameter 
as the mode of its probability distribution and the 1 a limits as the 
value at which the integral of the distribution equal s 34.1% from 
each side of the mode. The iGelman & Rubinl i ll 9921) statistic was 
calculated for each fitted parameter and we concluded that chain 
convergence was good. 

To test how the limb darkening coefficients affect the derived 
system parameters we experimented with fitting the limb darken- 
ing coefficients of the RISE light curves which have higher pre- 
cision. Our tests indicate that the quality of the RISE light curves 
is not enough to fi t for both limb darkening coefficients (see also 
iGibson et al]l2008l) . and hence, we choose to fit the linear LDC a 
which dominates. The quadratic LDC of the RISE light curves 6 
and both LDCs of the JGT light curve were kept fixed. Therefore, 
in this second MCMC procedure 21 parameters were fitted. 



4 RESULTS 

In Table [2] we present the system parameters of WASP- 13 and 
the lcr uncertainties derived from the MCMC analysis for both 
the fixed and fitted LDCs. The first set of parameters can be di- 
rectly derived from the light cur ve analysis. The parameters wer e 
derived using the equat i ons of ISeager & M allen-Orn elasl d2003l) : 
ISouthworth etai] d2007t) : iKippind <201Ch . Specifically, the Ti_4 
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Table 2. WASP- 1 3 system parameters. 



Parameter 



LDC fixed 



LDC fitted 



Linear LDC for RISE a 
Normalised separation a/R* 
Planet/star radius ratio R p /R t 
Transit duration Ti— 4 [hours] 
Transit duration Tri [hours] 
Orbital inclination / [degrees] 
Impact parameter b[R„] 
Stellar density p* [p ] 

Orbital semimajor axis a [AU] 
Stellar mass M* [M©] 
Stellar radius R, [R ] 

Planet mass M p [M Jup ] 
Planet radius R p [Rj up ] 
Planet density p p [pj] 



7.58 
0.09225332 



0.4255 

+0.15 
0.13 

+0.00080 
0.00081 
4 Q04O.O24 

C> - DU ' -0.015 

85.64 ± 0.24 
0.575 ± 0.020 
0.309 4 



-,+0.018 
-0.016 



0.337 ±0.033 

+0.14 
0.15 
+0.00083 
0.0010 

4.063°°^ 

3.548°g« 7 
85.19 ±0.26 
0.621 ±0.021 
a +0.019 
-0.015 



7.39 
0.09226 



0.288 1 



0537Q+ - 00077 

U.UOJ(9_ 00059 



1.085 ± 0.04 

r, + 0.031 

-0.041 



1.512" 1 



05362+ ' 00074 
U.UOdDZ_ 00085 

1.09 ±0.05 
1.559 ±0.041 



485+ - 042 
U.*OO_ 05 8 

i.ouj_q 062 

190+ - 028 

U.13U_ 014 



+0.044 

u -*' ' -0.049 

1 •3KQ+ ' 045 
1.009_ 055 

167+ ' 017 
u - lul -0.022 




6200 



6000 5800 5600 
Effective temperature (K) 



5400 



5200 



refer s to the total tra nsit duration and Tti is denned by equation 
15 of l Kipping! d2010h . 

The derived linear limb darkening coefficient, a = 0.34 ± 
0.03, is statistically significantly different from the theoretical 
value a — 0.4402. The chain convergence was good and the linear 
LDC parameter is well constrained by our observations. We adopt 
the fitted LDC solution because according to the F-test it is sig- 
nificantly (99.9%) better than the fixed LDC one. In section 5 we 
present a more in depth discussion about the discrepancy between 
the fitted and theoretical LDCs. 

All the parameters of both solutions are within ~ 3<r. We 
found that the inclination of the orbit is 85.2 ± 0.3 degrees and 
the derived stellar density is p* = 0.29 ± O.O2p . T his is signifi- 
cantly lower than previously assumed (p, = O.43p Skill en et al.l 
2009), implying an evolved star. To obtain the stellar and planetary 
radii we need to determine the stellar mass, as explained in the next 
section. 



4.1 Stellar mass and age 

Our high quality transit light curves allow an accurate estimate of 
the stellar density. This, combined with stellar models, can be used 
to constrain the stellar mass and age. We in terpolate the Yonsei- 
Yale evolution trac ks jDemarque et al.l |2004) for [Fe/H] = 
( Skill en et al.l l2009t) which are plotted in Figure [5] Also shown 
are the position of WASP- 13 with the effective temperature from 
ISkillen et alj d2009ri and our value of p„ . We derive a stellar mass 
of 1.09 ± O.O5M and a stellar age of 7.4 ± 0.4 Gyr. Both agree 
with the estimate in th e original dis c overy paper. Using the stel- 
lar reflex velocity ffom lSkillen et al.l {2009), Ki — 55.7 m/s, we 
estimate a planetary mass of 0.48 Mj up . 

In Figuref4]we show the evolutio nary tracks for a stellar mass 
of 1.0, 1.09 and 1.2 M adapted from lDemarque et al.l fe004h and 
we overplot the position of WASP- 13. These models suggest that 
WASP- 13 has evolved to the hydrogen- shell burning phase and is 
close to entering the sub-giant branch. However, the uncertainty 
in the WASP-13 metallicity ( [Fe/H] = ± 0.2) affects the de- 
ter mination of the evolutio nary status of WASP-13. According to 
fhe lDemarque et al.l |2004) models, if WASP-13 has a metallicity 



Figure 3. Isochrone models (sol id lines) fr om Demaraue et al. 1 2 0041) for 
WASP-13 using [Fe/H] =0 from I Skillen et alj 120091) and assuming solar 
composition Z=0.01 8 1 . The age in Gyr is marked in the left of the respective 
model. We also show the mass tracks (da shed lines) for 1.0, 1.05, 1.09 and 
1.2M . We overplot the T e/ f = 5826 ISkillen et al]|2009h and our fitted 
(p./po)" 1 / 3 for WASP-13. 




6000 5500 5000 

Effective temperature (K) 



Figure 4. Evolutionary mass tracks from Demarque et al. (2004) for the 
same stellar parameters as Fig.fJJfor stellar masses of 1.2, 1.09 and 1.0 M 
from left to right. For the 1.09 M evolutionary mass track we also show 
the 5, 6,7 and 8 Gyr points. 



of +0.2 then it could have a mass of 1.2 M which is above the 
critical mass for the stellar cores to become convective and where 
details on the treatment of convective overshooting became impor- 
tant. In this case, the evolutionary status would be more uncertain 
but probably WASP-13 would be in the contraction phase before 
the hydrogen- shell burning phase starts. 

We obtain a larger stellar and planetary radius (R, = 1.56 ± 
0. 04 R and K v = 1.39 ± 0.05 Rj up ) than previously reported 
in ISkillen et al.l d2009h . Our derived parameters are consistent with 
the "b=0.46" adopted solution presented in the discovery paper, al- 
though our precision is almost 10 times better. 
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Table 3. Time residuals from the linear ephemeris. 



Date 


Epoch 


Time residuals (sec) 


Uncertainty (sec) 


2008-02-16 


-244 


235 


82 


2009-01-29 


-164 


83 


42 


2009-05-05 


-142 


-211 


51 


2010-02-03 


-79 


-92 


39 


2011-01-13 





114 


44 



4.2 Transit times 

The estimated transit times were used to update the linear 
ephemeris, 



T(HJD) = T + EP. 



(1) 



We found P = 4.353011 ± 0.000013 days and T = 
2455575.5136 ± 0.0016, which was set to the mid-transit time 
of the 2011 light curve. The updated period was used in the final 
MCMC procedures. 

Time residuals from the linear ephemeris are given in Ta- 
ble [3] The partial transits show a larger deviation from the lin- 
ear ephemeris whi ch was also noted for the case of TrES-3b 
(Gibson et al. 2009). These authors concluded that the uncertain- 
ties in the transit times of partial transits are underestimated, prob- 
ably due to the lack of out-of-transit data needed for normalisa- 
tion. Therefore, although the time residual of the May 2009 transit 
is significant at 4 sigma, this could be due to normalisation errors 
rather than an intrinsic transit timing variation. For example, re- 
moving the out-of-transit data of the May 2009 light curve results 
in a difference in the estimated mid-transit time of 155 s = 3.5cr. 
This suggests the presence of additional red noise in the light curve 
that affects the transit times which should be regarded with cau- 
tion. Hence, we conclude that the time residuals of WASP- 13b are 
probably consistent with a linear ephemeris, but encourage further 
transit timing observations. 



5 LIMB DARKENING 

Our measured linear limb darkening parameter (a = 0.34 ± 0.03) 
is statistically different from the one predicted from stellar models 
(q = 0.4402). We c hoose the theoretical LDCs from stellar models 
of lHowarthl l l2011al) because they include the ' V+R' RISE filter. 

Disagreement between the empirical an d theoretical LDCs 
has been reporte d for other systems (e.g. South wortf] 1200 8t 
iBarros et al.l | 201 ll) including for higher quality light curves from 
HST dClaretll2009l). CoRoT dSindl2010l. Csizmadia. Sz. et al in 
prep) and Kepler dKipping & Bakosll201 ll) . This suggests that the 
discrepancy is not due to insufficient quality of our light curves. 
Moreover, Iciaretl d2009l) concluded that uncertainties in the stellar 
temperature and metallicity cannot explain the difference between 
the theoretical and empirical LDCs. Therefore, the causes of this 
discrepancy are not clea r. 

Recently, iHowarthl ( feOllbl) showed that LDCs derived from 
transit observations depend on the impact parameter and cannot be 
directly compared with the theoretical values. This is because for 
a higher impact parameter the transit shape is less sensitive to the 
global LD law. This can lead t o up to a 60% difference in the es- 
timated LDCs faowarthll2011bl) . Fortunately, this insensitivity also 
means that the effect on the transit parameters is relatively small. 



Using a synthetic photometrv/atmosphere-model. lHowarthl d20 lTbh 
predicts that this geometric effect leads to an underestimation of 
the linear LDC of the quadratic law which qualitatively explains 
our estimated linear LDC. If this was the only cause of the discrep- 
ancy we should adopt the theoretical LDC solution for the system 
parameters because it would remove the correlation between the 
LDCs and the impact parameter providing more accurate results. 

Using HST and Kepler data. lHowarthl (201 lbh concluded that 
accounting for this geometrical dependency leads to better agree- 
ment between the measured and predicted LDCs. However, this 
cannot reproduce all of the empirical data suggesting that there 
could be other causes of the discrepancy. Tentative explanations 
include the missing physics in the atmospheric models like, for ex- 
ample, the effects of faculae and granulation (Csizmadia, Sz. et 
al i n prep.) and the plane-p arallel approximation in stellar mod- 
els l lNeilson & Lesterll201 ll) . In this case, it would be advisable to 
adopt the fitted LDC solution for the system parameters because it 
is less dependent on the stellar models. 

Further investigation of very high precision transit light curves 
such as from the Kepler mission is needed to better understand the 
difference between the measured and theoretical LDCs. This will 
help to decide if we should fit for the LDCs in lower quality light 
curves. 



6 DISCUSSION AND CONCLUSION 

We have presented four new transit light curves of WASP- 13b in- 
cluding two full transits. To derive the parameters of the system, 
we used the same MCMC procedure as lBarros et alj ( 1201 11) by fit- 
ting these new light curves toge ther with a previous JGT transit 
of WASP- 13b dSkillen et al]|2009l) . We estimated the uncertainties 
accounting for formal errors and systematic noise. 

Our derived value of the linear LDC is statistically signifi- 
cantly different from the theoretical value predicted for stellar at- 
mospheres models. Although this does not significantly affect the 
derived system parameters it can be important for stellar atmo- 
sphere models. We adopt the fitted LDC solution because it is a 
statistical significantly better fit to the data. 

The impact parameter of WASP- 13b was poorly constrained 
in the discovery paper (0 < b < 0.46) though the higher value was 
adopted. We obtain b = 0.62 ± 0.02, which is close to the previous 
adopted value although slightly higher. This implies a larger stellar 
radius of R* = 1.56±0.04Rp) w hich combined with the evolution 
models of lDemarque et alJ d2004) points to an age of 7.4±0.4 Gyrs 
for the host star. It also suggests that WASP- 13 has evolved off 
the main-sequence and is now in the shell hydrogen burning phase 
which would imply that in less than 2 Gyr the star will engulf the 
planet's current orbit. The evolutionary status of WASP-13 favours 
a slightly higher mass for the host star of ~ I.IMq which also 
implies a slightly higher planetary mass of 0.48Mj up . 

We derived a larger planetary radius for WASP- 13b, (R p = 
1.39 ± 0.05Rj up ), than previously reported and hence a lower den- 
sity, (p p = 0.17 ± O.Olpj). Therefore, WAS P- 13b is the fifth low - 
est density e xoplanet known after K epler- 7b dLatham et al.ll2010h. 
WAS P- 17b dAnderson et"ai] |2010|). Tr ES-4b dMandushev et all 
120071) and CoRoT-5b dRauer et alJl2009T) and its properties are im- 
portant for irradiation models. We estimate an equilibrium temper- 
ature for WASP-13b of ~ 150 K. The hydrogen/helium coreless 
models of iFortnev et al.1 d2007l) predict a radius of ~ l.lRj up for 
WASP- 13b. This is in agreement with the radius derived from the 
polynomial fit of lLaughlin et al.l d201lh to the baseline models of 
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iBodenheimer et all d2003h . R p ~ 1.16. We derive a rad ius anomaly 
for WASP-13b of U = 0.23 jLaughlin et all 1201 lh . Therefore, 
WASP- 13b is inflated relative to the theoretical thermal irradiation 
models for gas giant planets. 

The circularization timescale for WASP- 1 3b is 0.0 1 1 -0. 1 1 Gyr 
for Q p = 10 5 — 10 6 . Therefore, we expect a circular orbit and it 
is unlikely that the inflated radius of WASP- 13 could be explained 
by tidal dissipation. Enhanced atmospheric opacitie s of 10 x solar, 
plus t he correction due to the transit radius effect dBurrows et al.l 
l2007t) . can explain a radius anomaly of ~ 0.18 for a 7.4 Gyr 
planet and hence might be able to explain the radius anomaly of 
WASP- 13b. The mod els of IBodenheimer et al J d2003h which in- 
clude kinetic heating (Guillot & Showman 2002) can also explain 
the radius of WASP - 13b. Finally, the ohmic dissipation model of 
Baty gin et al.l ilOl lh explains all known ex oplanet radii. 

According to the classification of iFortnev et al.l J2008), 
WASP-13b belongs to the highly irradiated "pM" class planets. 
These have considerable opacities due to molecular hazes of TiO 
and VO and are expected to have temperature inversions in their 
atmosphere. This can be tested by secondary eclipse observations. 
The large scale height of WASP- 13b, ~ 550 Km, and its bright host 
star, V = 10.4 , makes it also a good target for transmission spec- 
troscopy. Probing the atmospheric composition of this planet can 
shed light on its inflated radius. 
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